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For free fields, pair creation in expanding universes is associated with the building up of correla- 
tions that lead to non-separable states, i.e., quantum mechanically entangled ones. For dissipative 
fields, i.e., fields coupled to an environment, there is a competition between the squeezing of the 
state and the coupling to the external bath. We compute the final coherence level for dissipative 
fields which propagate in de Sitter space, and we characterize the domain in parameter space where 
the state remains non-separable. We then apply our analysis to (analogue) Hawking radiation by 
exploiting the close relationship between Lorentz violating theories propagating in de Sitter and 
black hole metrics. We establish the robustness of the spectrum and find that the entanglement 
among Hawking pairs is generally much stronger than that among pairs of quanta with opposite 
momenta. 



I. INTRODUCTION 

The propagation of quantum fields in expanding cos- 
mological backgrounds leads to the spontaneous creation 
of pairs of particles with opposite momenta [1]. For free 
fields, relativistic or dispersive, this pair creation (also 
called the dynamical Casimir effect in condensed mat- 
ter physics, see e.g., Refs. [2, 3]) is associated with the 
building up of non-local correlations that lead to quan- 
tum mechanically entangled states [4, 5] . To define these 
states without ambiguity, we shall use the notion of non- 
separability [6], see Appendix B. For dissipative fields, 
i.e., fields coupled to an environment, there is a compe- 
tition between the squeezing of the state which increases 
the strength of the correlations, and the coupling to the 
external bath which reduces it [7-9]. 

Our principal aim is to study this competition. We 
shall work both in time-dependent (cosmological) set- 
tings and with stationary metrics. For simplicity and def- 
initeness, we consider fields which propagate in de Sitter 
space, and which display dissipative effects above a cer- 
tain momentum threshold A. For these fields, the final 
coherence level is constant and well defined. We char- 
acterize the domain in parameter space where the final 
state is non-separable. The parameters are mass of the 
field, the temperature of the environment, and the ra- 
tio A/if, where H is the Hubble constant. Since the 
dissipative/dispersive effects we are considering are sup- 
pressed in the infrared, our models can be conceived as 
providing a phenomenological approach to theories of 
quantum gravity, such as Hofava-Lifshitz gravity [10], 
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where Lorentz invariance is violated at high energy. In 
these theories, dissipative effects will necessarily appear 
through radiative corrections [11]. We also recall that in 
condensed matter, the spectrum of quasi-particles often 
displays dissipation above a certain threshold. Hence, 
our model can also be viewed as a toolbox to compute 
the consequences of dissipation on pair production and 
parametric amplification found, e.g., in the superfluid of 
polaritons studied in Ref. [12]. 

The interest of working in de Sitter space is twofold. 
On the one hand, the analysis of the state can be done in 
terms of homogeneous modes and pair creation of quanta 
with opposite momentum. On the other hand, the state 
can also be analyzed in terms of stationary modes and 
thermal-like effects associated with the Gibbons-Hawking 
temperature [1]. ft is rather clear that the homogeneous 
representation in de Sitter can be conceived as an ap- 
proximation to e.g., slow roll inflation, see Refs. [13, 14]. 
What is less obvious is that de Sitter also provides a re- 
liable approximation to describe dissipative fields propa- 
gating in black hole metrics. Indeed, when the ultravio- 
let scale A is well separated from the surface gravity of 
the black hole, the dissipative aspects of typical Hawking 
quanta all occur in the near horizon region which can 
be mapped into a portion of de Sitter space (when the 
Hubble constant is matched to the surface gravity). As a 
result, the state evaluated in a black hole metric can be 
well approximated by the corresponding one evaluated in 
de Sitter. In this respect, the present paper follows up 
on our former work [15] where we studied this correspon- 
dence for dispersive fields. The reader unfamiliar with 
field propagation in de Sitter space will find in that work 
all necessary information. 

This paper is organized as follows. In Section II we 
present the action which engenders dissipative effects, 
and we discuss the residual symmetries found in de Sitter 
space when considering such theories. In Section III, ex- 
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ploiting the homogeneity of de Sitter, we compute the 
spectral properties and the correlations of pairs with op- 
posite momenta. In Section IV, exploiting the station- 
arity, we compute the deviations with respect to the 
Gibbons-Hawking temperature. We apply our model to 
black holes in Section V, and we conclude in Section VI. 
We work in units where h = c = 1. 



II. DISSIPATIVE AND DISPERSIVE FIELDS 
A. Covariant settings 

We study a 1 + 1 dimensional scalar field cp which 
has a standard relativistic behavior at low energy but 
which displays dispersion and dissipation at high en- 
ergy, thereby violating (local) Lorentz invariance. While 
high energy dispersion is rather easily introduced and 
has been studied in many papers both in cosmological 
settings [14, 16-18] and black hole metrics [19, 20], see 
Ref. [21] for a review, dissipation has received compara- 
tively much less attention and is technically more difficult 
to handle when preserving unitarity and general covari- 
ance. To do so in simple terms, following Refs. [13, 22], 
we introduce dissipation by coupling <p to some environ- 
mental degrees of freedom ip, and the action of the entire 
system Stot = S^ + S.0 + S- m t is taken quadratic in (p,ip, 
as in models of quantum Brownian motion [23]. 

In the present work, we consider phenomenological dis- 
persion relations which contain both dispersive and dis- 
sipative effects. These relations can be parametrized by 
two real functions T, / as 



n 2 + 2iTfl = 



P 2 + f = F 2 



(1) 



where T(P 2 ) > is the damping rate, and f{P 2 ) de- 
scribes dispersive effects. To recover a relativistic behav- 
ior in the infrared, a typical behavior would be T ~ P 2 
and / ~ P 4 for P 2 — > 0. In Eq. (1), Q and P 2 are, respec- 
tively, the proper frequency and the proper momentum 
squared as measured in the "preferred" frame [24], i.e., 
the frame used to implement the dispersion relation. In 
condensed matter systems, it is provided by the medium. 
Instead, in the phenomenological approach to Lorentz vi- 
olating effects we are pursuing, it should be given from 
the outset, either as a dynamical field endowed with an 
action [10, 25], or as a background field (as we shall do). 
To describe it in covariant terms, following Ref. [26], we 
introduce both the unit timelike vector field u which de- 
scribes the flow of preferred observers, and the unit space- 
like vector field s which is orthogonal to u. In terms of 
these, one has = u^p^ and P 2 = (s^p^) 2 where p^ is 
the momentum of the particle in an arbitrary coordinate 
system. In two dimensions, the metric g M1/ can be written 
as g M „ = —u^u v + s^Sv which expresses that u and s are 
orthonormal vectors. 

We now consider a unitary model which implements 
Eq. (1). This model is not unique but can be considered 



as the simplest one, as shall be made clear below. In 
covariant terms, the total action S'tot = S<f, + + S- ln t is 

S'tot =\ Jd 2 [-g^V^V^-mV- 

+ \fd 2 jdq [{V u %f - (nAqf ^ 2 ] (2) 



7(V S )0 V u dqip q 



where d 2 = d 2 x^J —g{x) is the covariant measure. In 
the first line, S^ is the standard action of a massive 
scalar field, apart from the last term which introduces the 
high frequency dispersion described by f(P 2 )- The self- 
adjoint operator which implements P 2 is — = VjV s , 
where V s = s M V^ is an anti-self adjoint operator (when 
u is a freely falling frame), VJ = — V^s^ its adjoint, and 
V K the covariant derivative. 

The second action, that of the ip field, contains the 
extra dimensionless parameter q, which can be consid- 
ered as a wave number in some extra dimension. Its 
role is to guarantee that the environment degrees of free- 
dom are dense, something necessary to engender dis- 
sipative effects when coupling ip to <p [22, 23]. The 
role of the frequency A is to set the ultraviolet scale 
where dissipative effects become important. The kinetic 
term of ip is governed by the anti-self adjoint operator 
V u = — (u^Vp + V M w M )/2 which implements fl = u^p^. 
We notice that there is no spatial derivative acting on 
ip. This means that the quanta of ip are at rest in the 
preferred frame. This restriction can easily be removed 
by adding the term c^CVgip) 2 which associates to the 
group velocity of the low q quanta. Since the inclusion 
of this term would lead to more complicated equations, 
without adding more physics, we work with = 0. 

The interaction between the two fields is given by the 
third action. The strength and the momentum depen- 
dence of the coupling is governed by the function j(P) 
which has the dimension of a momentum. Its role is to 
engender the decay rate T entering Eq. (1). The last 
two actions possess peculiar properties which have been 
adopted to obtain simple equations of motion. These are 



[V/u"V, + F 2 (~V 2 S )] 4> = 7 (V t s )V„ J dq^ q , (3a) 
[Vl + (nAq) 2 ] ^ = -V„ 7 (V s )0 . (3b) 
The solution to the second equation is 



^ ? (x') = ^°(x') - / d 2 G g (x', x)V u7 (V 5 Mx) , (4) 

where ip® is a homogeneous solution, and where the 
driven solution is governed by G q (x, x'), the retarded 
Green function of ip q . When injecting tp q in the r.h.s. of 
the first equation, one obtains the equation of (p driven by 
ip q . The general solution can be written as (p = (p dcc +<p dT , 
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where the decaying part is a homogeneous solution, and 
where the driven part is given by 



#V) =/ d 2 G ret (x',x) 7 (Vt)V» / d#°(x) 



(5) 



In a general Gaussian (f> — ip model, the retarded Green 
function G rct would obey a nonlocal equation, i.e., an 
integro-differential equation. We have adjusted the prop- 
erties of S^p and S- ln t precisely to avoid this. Two proper- 
ties are essential. Firstly, at fixed q and along the orbits 
of u, Eq. (3b) reduces to that of a driven harmonic os- 
cillator. This can be seen by introducing the coordinates 
(t, z) defined by u^d^ = —d T \ z where z is a spatial coor- 
dinate which labels the orbits of u. Then, V u applied on 
scalars is 



(6) 



where a(r, z) = eJ dT ' e ( T < z ) ( and where 9 = — V M w M is 
the expansion of u. Hence the rescaled field 

*g {r,z) = y/a{T, z) ip q (t, z) (7) 

obeys the equation of an oscillator of constant frequency 
7rA|g|. Secondly, when summed over q, the retarded 
Green function of ip obeys [22] 



/CO 
dqG q (x,x r ) 
-OO 



g(x-xQ 
A 



(8) 



where <5 2 (x — x') is the covariant Dirac delta, i.e., 
/ d 2 f{x)8 2 (x - x') = /(x'). Eq. (8) guarantees that 
the differential operator encoding dissipation is local. 
Namely, when inserting i[i q of Eq. (4) in Eq. (3a), one 
finds 



□di S8 = 7(V|)V„ J dq^° q , 



(9) 



with the local differential operator 



□ 



diss ■ 



7(Vt)V„ 7 (V.) 



A 



(10) 



One can now verify that the WKB solutions of 
Ddiss 4> = are governed by a Hamilton- Jacobi action 
which obeys the dispersion relation of Eq. (1) with 



r = M 2 /2A, 



(11) 



see Appendix C for more details. The reader can also 
verify that any modification of the actions and Si n t 
leads to replacement of Ddiss by a non-local operator. 
When considered in homogeneous and static situations, 
this is not problematic because one can work with Fourier 
modes in both space and time. However when considered 
in non-homogeneous and/or non-static backgrounds, it 
becomes hopeless to solve such an equation by analytical 
methods. 



In our model, the retarded Green function thus obeys 

□ diss G rc t (X,X') =5 2 (X-X') , (12) 

and vanishes when x is in the past of x', where the past is 
defined with respect to the foliation introduced by the u 
field. When canonically quantizing (j) and tp, since our ac- 
tion is Gaussian, the commutator G c (x, x') = [0(x), cj)(x')] 
is independent of /5 to t, the state of the entire system. 
Moreover, it is related to G re t in the usual way 



-iG c (x, x') — G ret (x, x') — G ret (x',x) 



(13) 



In this paper we only consider Gaussian states. This 
implies [27, 28] that the density matrix /5 to t, and all 
observables, are completely determined by the anti- 
commutator of d>, 



G ac (x,x')=Tr Ptot{<Kx),<£(x')} , 



(14) 



that of and the mixed one containing (j) and tp. De- 
composing the field operator <fi = dcc + </> dr , G ac splits 
into three terms. The first one involves only (j> dcc , the 
second contains both (j> doc and (fr T , and the last only dr . 
When assuming that the initial conditions are imposed 
in the remote past, because of dissipation, only the last 
one is relevant. Using Eq. (5), it is given by 

G d J(x,x')= JJd 2 1 d 2 2 G Ict (x,x 1 )G rct (x',x 2 )N(x 1 ,x 2 ) , 

(15) 

where the noise kernel is 

iV(x,x')=7(V t s )V tl7 (V?)VL 



dqdq'Tr(p to S°(x),4°Ax')} 



(16) 



In Sections III and IV, we compute Gf J and extract from 
it pair creation probabilities and Hawking like effects tak- 
ing place in de Sitter space. 



B. Affine group in de Sitter space 

The two dimensional de Sitter space possesses three 
Killing vector fields that generate the algebra of the Lie 
group SO(l, 2). Imposing that the action is invariant un- 
der the full group precludes ultraviolet dispersive and dis- 
sipative effects such as those of Eq. (1), see Appendix A 
in Ref. [15] for the proof. Since we want to work with 
Eq. (1), we must break (at least) one of these symme- 
tries. As in Rcfs. [15, 29], we preserve the invariance 
under a two dimensional sub-group which corresponds to 
the affine group. Its algebra is generated by the Killing 
fields K z and K t . Using the cosmological coordinates t, z 
of the Poincare patch 



ds 2 



-dt 2 



„2Ht 



dz 2 



1 



H 2 ^ 



[-drf + dz 2 



(17) 
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K z — d z generates translations in z and expresses the 
homogeneity of the sections t = est., whereas K t — dt — 
Hzd z expresses the stationarity of de Sitter. Using X — 
e z, this symmetry becomes manifest, 

ds 2 = -dt 2 + (dX - HXdtf . (18) 

Considering the two Killing fields K z and K t , there is 
only one unit time-like freely falling field which commutes 
with both of them. We call it u s , and we call sg the 
spatial unit orthogonal vector u s ■ Sff = 0. Then the 
coordinates t, X are both invariantly defined in terms of 
u s , s ff by dt = uf t dx^, d x \t = s^d^. 

Imposing that the action of Eq. (9) be invariant under 
the affine group, requires that the preferred fields u and 
s commute with K t and K z . This fixes u and s up to a 
boost, see Appendix A. For simplicity, in what follows, 
we work with u — u s . In this case, the preferred frame 
coincides with the cosmological one, and the orbits of u 
are z = est. 

We also impose that the states ptot are invariant under 
the affine group. This is analogous to the restriction to 
the so-called a— vacua which are invariant under the full 
de Sitter group [30, 31]. This means that G ac , G lct and N 
of Eq. (15) will be invariant under both K t and K z . How- 
ever, because the commutator [K z ,K t ] = —HK Z does 
not vanish, one cannot simultaneously diagonalize K z 
and K t . This leads to two different ways to express the 
2-point functions, either at fixed wavenumber k = — id z \ t , 
or at fixed frequency lo = id t \x- Explicitly, one has 

GL y (v, rf) = J dAz e~* kAz G any (Az, V , rf) , (19a) 
GZ ny (X,X>) = J dAte^ At G any (At,X,X') , (19b) 

where k = |k|, and where the "any" subscript indi- 
cates that these Fourier transforms apply to any 2-point 
function which is invariant under the affine group. (In 
Eq. (19a), G k only depends on k because we impose 
isotropy.) 

What is specific to this group is that the two symme- 
tries combine in a non trivial way, and imply that 2-point 
functions only depend on two quantities, and not three, 
as it is generally the case in homogeneous or stationary 
metrics. In the homogeneous representation, it implies 
that the product k G any (77, rf) only depends on the phys- 
ical momenta P = —Hkrj, P' = —Hkrj'. Hence, in what 
follows, we work in the P-representation with 

G any (P,P') = ^Gi iy (r ) ,r,') . (20) 

To reach this representation when starting from the sta- 
tionary G^ ay (X, X') is more involved, and is explained 
in Appendix A. 



III. HOMOGENEOUS PICTURE 
A. Dissipation and non-separability 

In this section, we decompose the fields in Fourier 
modes of fixed k. This representation is suitable for 
studying the cosmological pair-creation effects induced 
by the expanding scale factor a(t) — e Ht = —\jHr\. 

To express the outcome of dissipation in standard 
terms, we exploit the fact that Lorentz invariance is re- 
covered in the infrared, for momenta P — ke~ Ht -^i A. In 
this limit, since T and / of Eq. (1) are negligible, the k- 
components of 4> decouple from tp, and obey a relativistic 
wave equation. Hence, the k-component of the (driven) 
field operator of Eq. (5) can be decomposed in the out 
basis as 

Mt) ~ a k ^(t)+a f _ k ^(t) , (21) 
t— ^00 

where the out modes obey the scalar wave equation and 
satisfy the standard positive frequency condition at late 
time. This means that the (reduced) state of (f> (ob- 
tained by tracing over ip) can be asymptotically described 
in terms of conventional excitations with respect to the 
asymptotic ouf-vacuum. 

The out operators a,k, a^ obey the standard commuta- 
tion rule [at, a^,] = S(k — k'). For notational simplicity, 
we omit the S(k — k') when writing two-point functions 
because it is common to all of them since we only consider 

homogeneous states. For instance, Tr^ptot {^k, 4>k'}j = 

S(k - k') x G k ac . Using Eq. (21), the coefficient of the 
(5-function is 

G ~(*.*)Loo = 2 ^ 2nk + \Mt)\ 2 + 4Mcwl(t)) , 

(22) 

where 

n k = Tr^ptotaj^ak) = Tr^ptot a^ k a_ k ) , (23a) 

c fe = Tr(p tot a_ k ak) ■ (23b) 

The mean number of asymptotic outgoing particles is 
n& > 0, whereas the complex number Cfc characterizes 
the strength of the correlations between particles of oppo- 
site wavenumber. The relative magnitude of this number 
leads to the notion of non-separability. 

To explain this, we recall that the correlations weighted 
by Cfc obey the following Cauchy- Schwartz inequality, 

|cfc| 2 <nfc(n fe + l) , (24) 

see Appendix B for more details. To characterize the 
level of coherence, we shall use the parameter of Ref. [32] 

5k =n k + l- \c k \ 2 /n k , (25) 

which belongs to the interval [0,nfc + 1]. When 8 k = 0, 
one has a maximally entangled squeezed state with zero 
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entropy, and when 8k = Uk + 1 one has an incoherent 
thermal state of maximum entropy. For homogeneous 
Gaussian states, one also verifies that the entropy is 
monotonically growing with 8k- 

The important and non-trivial fact is that 8k = 1 di- 
vides states that are quantum mechanically entangled 
from states that only possess classical correlations. To 
show this we recall the notion of separability. A two- 
mode state is called separable when it can be written 
as a weighted sum of products of two one-mode states, 
where all weights are positive and can thus be interpreted 
as probabilities. In this case, the strength of the correla- 
tions is more restricted than Eq. (24). Indeed, one finds 
l c fc| 2 < n\, see Appendix B. As a consequence, whenever 



n\ < \c k \ 2 < n k (n k + 1) , 



(26) 



a homogeneous state is non- separable, i.e., so entangled 
that it cannot be represented as a classically correlated 
state characterized by probabilities. In terms of 8k this 
criterion is simply given by 8 k < 1. 



B. Invariant states and P-representation 

Since the states we consider are invariant under the 
afHne group, n k and Cfc are necessarily independent of 
k. We shall nevertheless keep the label k to remind the 
reader that we work at fixed k and not at fixed u as in 
the next section. Because of the affine group, 



<p(P) = sfkfE x <p k (t) , 



(27) 



only depends on P, where ipk{t) is the (positive unit 
norm) out mode of Eq. (21). The norm of the mode 
ip is fixed by the Wronskian 



W{ip) = 2H 2 Im(cp*d P (p) = 1 



(28) 



Using such ip and Eq. (20), Eq. (22) and Eq. (15) can be 
written as 

GUP P)\p^o = 2 I 2 ™* + 1] \<f{P)\ 2 + {cW 2 (P)) , 

(29a) 

D dP 1 dP 2 _ 

p2 p2 
'0 r l r 2 



p2 p2 G rct (p Pi)G iet (P, P 2 )N(P 1 ,P 2 ) 



(29b) 

In the second line, the noise kernel of Eq. (16), which 
is also invariant under the affine group for the set of 
states we are considering, has been written in the P- 
representation using Eq. (20). To extract rik and c k from 
the above equations, we need to compute G ret and N. 
Using Eq. (20), Eq. (12) reads 



H 2 2 ^ 



aVp 



Hd 



7 VP P 2 



G Iet (pP')=8(P-P'). 

(30) 



The unique (retarded) solution can be expressed as 

G rct (P, P') = 26{P' - P)Im (<p P <p*p, ) e- x ? , (31) 
with the optical depth [13] 



(32) 



Its role is to limit the integrals over Pi and P 2 in 
Eq. (29b) to low values so that Xq < 1. All informa- 
tion about the state for higher values of P is erased by 
dissipation. In Eq. (31) we have introduced 



(p P = e x ° ip(P) , 



(33) 



where (p is a homogeneous damped solution of Eq. (30) . 
By construction, Cpp obeys the reversible (damping free) 
equation 1 



[H 2 P 2 d 2 P + F 2 - T 2 ] (p P = 



(34) 



and is normalized by Eq. (28). Moreover, we impose that 
it obeys the out positive frequency condition, meaning 
that in the limit P — > 0, it asymptotes to the out mode 
ip of Eq. (27). Hence, comparing Eq. (22) with Eqs. (29b) 
and (31), we find 

(35a) 
(35b) 



' dPl dP 2 ; ' I 

p2 p2 
r l r 2 



(pp^e-^ 1 -^ 2 N(Pt,P 2 ) 



These central equations establish how the environment 
noise kernel N fixes the late time mean occupation num- 
ber and the strength of the correlations. 

We now compute N. When u is freely falling, the 
rescaled field ^1°. of Eq. (7) is a dense set of independent 
harmonic oscillators of frequency fl q = 7rA|g|, one at each 
z. This simple setting arises because we set = in 
the action for i/j, see the discussion after Eq. (2). It also 
means that the positive frequency mode functions are 
the standard e~ tfi «* / \/2Q q , and that the state of these 
oscillators can easily be described in a Fock-space decom- 
position. 

In what follows we shall work with homogeneous ther- 
mal states of temperature T^. This means that the ex- 



1 For high values of P, the effective dispersion relation is super- 
luminal if dp(f — T 2 ) > 0, and subluminal if this quantity is 
negative. The critical case, / — T 2 =0, gives rise to a relativistic 
dispersion. In the case where F 2 — T 2 becomes negative, the 
mode enters an overdamped regime, see Ref. [13]. To avoid the 
complications this entails, we will only consider / — T 2 > 0. 
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pectation value of the anti-commutator of ip® is given by 

5(z-z') 



Tr(p tot {<(x),^(x')}) = 

v ' V a w a ( t ) 



x coth 



£lq COS (fl q At) 



2Tii, 



9.,, 



(36) 



The factor coth(f2 g /2T^) = 2n* + 1 is the standard 
bosonic thermal distribution. The prefactor S(z — 
z') I yj a{t)a(t') comes from the facts that *° of Eq. (7) 
is a dense set of independent oscillators, and that a(r, z) 
reduces here to the scale factor a(t). To get TV of Eq. (16) 
one should differentiate the above and integrate over q. 
The integration gives a distribution which should be un- 
derstood as Cauchy principal value: 



dqdq'V u V u ,Tr(p to S M^')}) = ~ ^AJjl 
v ' \j a(t)a(t' 



x i^ p - v - C0th(7r7 > Ai) 



) 

(37) 



To be able to re-express Eq. (37) in the P-representation, 
it is necessary to verify that it is invariant under the affine 
group. This is easily done using notations of the Ap- 
pendix A. One verifies that the first factor simply equals 
S(A2), whereas the second line is only a function of Ai. 
Taking into account the derivatives of Eq. (16), in the P- 
representation, the noise kernel at temperature Tw, reads 



d , (nT* , P' 

T P.V. coth — ^ln — 
1 H P 



91n^ 



(38) 



The symbol P.V. indicates that when evaluated in the 
integrals of Eq. (35), the non-singular part should be ex- 
tracted using a Cauchy principal value prescription on 
ln(P'/P) = H(t-t'). 

In the high temperature limit, the double integrals of 
Eq. (35) can be evaluated analytically because N ef- 
fectively acts as a Dirac delta function. Instead, when 
working with an environment in its ground state, or 
at low temperature T.^, we are not aware of analytical 
techniques to evaluate these integrals. Hence, to study 
the impact of dissipation on coherence in (near) vacuum 
states, we shall numerically integrate Eqs. (35). 



C. Numerical Results 

In the forthcoming numerical computations, for sim- 
plicity, we work with 



P 4 



,P 2 




n k for g 2 = 2 
■■ 5 k for g 2 = 2 
n k for g 2 = 



i±JL_ 
10 1 



iiiL 
10 2 



1 A 2 ' r 9 ~ 2A ' 



(39) 



FIG. 1. Numerical values for n k and 5 k for a massless field 
with critical damping g — g CT1 t and quartic superluminal dis- 
persion at the energy scale A = HX. For comparison, we have 
represented by a dotted line the n k of the quartic dispersive 
field (in which case 8 k = identically). Surprisingly, the state 
is non-separable, S k < 1, for all values of A. Moreover, 5k de- 
creases when dissipation increases. 



which contain the same ultraviolet momentum scale A. 
The dimensionless coupling g 2 controls the relative im- 
portance of dispersive and dissipative effects. In the limit 
g 2 — > 0, we get the quartic superluminal dispersion stud- 
ied in Refs. [14, 15]. The critical coupling g 2 rit = 2, 
greatly simplifies the calculations, since / — T = guar- 
antees that (pp obeys a relativistic equation, see Eq. (34). 

Using a numerically stable procedure for the Cauchy 
principal values like in Ref. [13], we compute and Cfe of 
Eq. (35) in the parameter space A, g 2 , m 2 , and T^. Since 
all physical effects only depend on dimensionless ratios, 
we present the numerical results in terms of /j. — m/H, 
A = A/H, and ?? = T^/H. 



1. Massless critical case 

We begin with the massless case (fi 2 — 0) and with 
g = g CI it- Then Eq. (34) is particularly simple since the 
rescaled mode Hp of Eq. (33) reduces for all P to the 

H. In this we recover the con- 



cwi-mode ipp — e lP / H / 
formal invariance of the massless field in two dimensions. 
There usually would be no particle production when it 
propagates in de Sitter space, however, the conformal in- 
variance being broken by dissipation, pair-creation will 
take place. 

In Fig. 1 we present rik and 8k when the environment 
is in its ground state (T^ = 0). For comparison, we also 
show rifc for quartic dispersion (g 2 — 0) which can be 
computed analytically in the Bunch-Davies vacuum [14]. 
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FIG. 2. Contour plots of lnn^ and \ri8k for a massless field with critical coupling g = g cr it in the parameter space (A = 
A./H, ■& = T-^/H) . At low temperatures, for i? = T^/H < f/10, and Sk barely depend on On the contrary, for high 
temperatures, $ > 1, n*, scales as oc $A~ 1//2 whereas 5fe scales as Sk oc tfA" 1 / 2 for A > 1, and S k oc i9A 1/2 for A < 1. The 
hatched region indicates the numerical uncertainty about the threshold value 5k = 1 found when nt « 1. 



For A — > co the number of particles goes to zero as 1/A, 
as is expected since conformal invariance is restored in 
this limit. Despite dissipation, we find that Sk < 1 for all 
values of A. This indicates that the state is always non- 
separable in the two-mode fc-basis. In addition, contrary 
to what might have been expected, the two-mode entan- 
glement is stronger for smaller values of A, i.e., stronger 
dissipative effects. The reason for this has to be found 
in the fact that A also sets the scale where conformal 
invariance is broken. 

Let us now turn to the effect of the environment tem- 
perature T^. Before presenting the results, we recall that 
for relativistic and dispersive fields, the vacuum state of 
zero temperature is the only stationary state which is 
Hadamard [15]. Hence the temperature is fixed to zero. 
This is not the case in our model: any temperature is 
acceptable, because the wave vector of the ip modes is not 
redshifted by the expansion (since we work with = 0) . 
As a result, the state of <j) can be heated up in a station- 
ary manner by the coupling to tp. Fig. 2 shows contour 
plots of rife and 8k for a massless field with Eq. (39), again 

for 9 = ffcrit- 

First, in the limit A — » oo, we observe that rik — > 
irrespectively of the value of T$. This establishes that 
there is a robustness of the relativistic result in the limit 
A — > oo which generalizes that found for dispersive fields, 
see e.g., Ref. [14]. Moreover, in the high temperature 
limit, Eqs.(35) can be evaluated analytically to give 



n k + \ 



6 k 



s/X 
s/X 



j 2 A 



(40a) 
(40b) 



where erfi is the imaginary error function. We com- 
pared the corresponding contours with the numerical 
ones shown in Fig. 2 and found that they are practically 
indistinguishable for i) > 10. 

Second, when considering the effect of the temperature 
of the environment, we observe two regimes. At low tem- 
perature (i? <C 1), n-k and 8k only depend on A and are 
basically given by the zero temperature limit shown in 
Fig. 1. However, at large temperature (•& ^> 1), they de- 
pend on A and $ according to Eqs. (40). As expected, the 
strongest signatures of quantum entanglement, 8k <C 1, 
are found in the region where the breaking of conformal 
invariance is large (and hence pair-creation is active) and 
when the environment temperature is small, so that the 
spontaneous pair-creation events are not negligible with 
respect to thermally induced events. On the other hand, 
when the temperature is large, the final state is separa- 
ble since Sk ^S> 1. In Fig. 2 (right panel) we see that 
the threshold case 8k = 1 is approximatively given by 
■& - X- 1 ' 2 for A < 1. The hatched region for A > 10 
represents the numerical uncertainty in the region where 
rifc is much smaller than 1. 



2. Massive fields 

We note that the massless case /i 2 = is an isolated 
point in the mass spectrum: a well-defined notion of out- 
quanta requires either fi 2 = or /i 2 > 1/4. In the latter 
case, the asymptotic owt-modes with positive frequency 
(see, e.g., Appendix B of Ref. [14]) are given by 



2 sinh n(i H 



(41) 




A 



FIG. 3. Contour plots of lnn^ and ln5k for a massive field (/i 2 
(A = A/H, -ff = T^/H). As in Fig. 2, for low temperature ■& < 
# > 1 and A 3> 1, rife and <5fc scale both as i9A~ 1//2 . 

where /2 = \j [i 2 — 1 / 4 and J denotes the Bessel function 
of the first kind. 

Fig. 3 shows the contour plots of rik and 5k for a mas- 
sive field with /1 2 = 5/4 and g = <? C rit, in the same pa- 
rameter space (A, i9) as in Fig. 2. The case of a Lorentz- 
invariant field in the Bunch-Davies state is recovered in 
the limit A — > 00, i9 — > 0. Now conformal invariance is al- 
ready broken by the mass term and therefore rik remains 
non-zero in this limit. 

At zero temperature, the strongest entanglement (low- 
est 5k) is found at large values of A, i.e., weak dissipation. 
This was expected, since dissipation reduces the strength 
of correlations. However, as in the massless case, the 
threshold of separability 5k — 1 is not crossed. 

When increasing the environment temperature T^, we 
see that the strength of correlation is reduced, and sep- 
arable states are found. The non-separability criterion 
5k < 1 is therefore only met either when is smaller 
than the Gibbons-Hawking temperature Tqh = H/2tt, 
or when the coupling to the environment is sufficiently 
weak. Notice also that the behavior at high temperature 
can again be obtained analytically, the integrals over the 
Bessel functions becoming hypergeometric functions. 

3. Role of g in the underdamped regime 

It is also interesting to consider the role of the coupling 
g, see Eq. (39). As g 2 approaches zero, the dissipative 
scale 2A/g 2 is moved deeper into the UV with respect to 
the dispersive scale which is fixed by A. In the limit g 2 — > 
0, the field becomes purely dispersive and rife, 5k can be 
computed analytically [14] in the Bunch-Davies vacuum. 
For g 2 < 2 the mode is underdamped. In this case, the 
solutions to Eq. (34) which correspond to asymptotic out- 




let* 10' 10° 10 1 10 2 

A 



= 5/4) with critical coupling g 2 rit = 2 in the parameter space 
.1, rifc and 5k are independent of the temperature. Instead for 



modes of positive frequency are given by, see Appendix 
BofRef. [14], 

g-TTA/4 R / p2 \ 

VP V A* 2 ' 2 V 2XH 2 J 

where A — g 4 and M is a Whittaker function 

defined in Ref. [33]. 

Fig. 4 shows contour plots of 5nfe/n° = (rife — n>k)/ n k 
(where n° is the number of particles without dispersion 
and dissipation) and 5k for a massive field in the under- 
damped regime. Here, we set = 0, and plot the results 
in the parameter space spanned by the two (dimension- 
less) ultraviolet scales: A which characterizes dispersion, 
and 2A/g 2 which is the UV scale of dissipation. The lat- 
ter is larger than the former in the underdamped regime. 
The grey areas therefore correspond to the overdamped 
regime which we did not study. 

In the weak dispersive/dissipative regime A > 10, it 
is evident that 5nk and 5k are both dominated by dissi- 
pative effects. For the latter, this is because dispersion 
alone does not lead to decoherence. For the deviation 
5nk , this follows from the fact that dispersion gives an ex- 
ponentially small correction to the pair creation process 
(see Ref. [14]), while the corrections due to dissipation 
are only algebraically small. As a result, the hierarchy of 
scales does not directly fix the importance of the respec- 
tive effects. 

On the other hand, when dispersion is strong (A < 1) 
the pair creation process is basically governed by disper- 
sive effects. The correction to the particle number due 
to dissipation is very small (compared to the dispersive 
correction) . One can also observe that the degree of two- 
mode entanglement is then basically governed by the sep- 
aration between the two scales g 2 , i.e., 5k is determined 
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FIG. 4. Contour plots of 5rife/n° and <5fc for a massive field (/1 2 = 5/4) in the underdamped regime g 2 < g^a- The environment 
is in its ground state ($ = 0) and the two axes are the dispersive scale A and the dissipative one 2X/g 2 > A. 



by the strength of dissipation at the dispersive threshold, 

(r/P) |p=a. 



IV. STATIONARY PICTURE 

In the absence of dispersion/dissipation, it is well 
known that the Bunch-Davies vacuum is a thermal 
(KMS) state at the Gibbons-Hawking temperature 
Tgh = H/2ir [1]. It is also known that this is the temper- 
ature seen by any inertial particle detector, and that this 
is closely related to the Unruh effect found in Minkowski 
space, and to the Hawking radiation emitted by black 
holes [34]. In the presence of dissipation, while the sta- 
tionarity of the state of (f> is exactly preserved when the 
state of the environment is invariant under the affine 
group, the thermality of the state is not exactly pre- 
served. This loss of thermality, which generalizes what 
was found for dispersive fields [15], questions the status 
of black hole thermodynamics when Lorentz invariancc 
is violated [35-37]. 



A. Loss of thermality 

To probe the stationary properties of the state, we con- 
sider the transition rates of particle detectors at rest with 
respect to the orbits of K t . This means that the detec- 
tor is located at fixed H\X\ < 1 in the coordinates of 
Eq. (18). In this case, the 2-point functions only depend 
on t — t' and can be analyzed at fixed u> = idt\x, see 
Eq. (19b). (The above restriction on X simply expresses 
that the trajectory be time-like.) 

The transition rates are, up to an overall constant, 
given by Fourier transforms of the Wightman function 
Gw [34]. The rates then determine n u (X) : the mean 



number of particles of frequency ui > seen by a detector 
located at X, through 



n»(X) = Gfy(X,X) 
UX) + l G W "(X,X) 



(43) 



To study the deviations with respect to the Gibbons- 
Hawking temperature Tgh = H/2n, we introduce the 
temperature function T U (X) defined by 



n u (X) + l 



-w/T u {X) 



(44) 



It gives the effective temperature seen by the detector, 
and reduces to the standard notion when it is indepen- 
dent of a;. In the absence of dispersion and dissipation, 
T U (X) = Tgh for all values of u>, which means that the 
Tolman law is satisfied [15]. 

In the following numerical computations, for simplicity, 
we work at X = with an inertial detector, with g — 
5crit: m — 0, and Tp = 0. Since the calculation of the 
commutator of (j> is much faster and more reliable than 
that of the anti-commutator, instead of using Eq. (43), 
shall be computed with 



>(X) 



G 



W (X, X) 
G«(X,X) 



(45) 



The denominator is expressed using Eq. (13). The nu- 
merator is obtained from Eqs. (69) and (38) with —>■ 0. 
In addition, the principal value is replaced by a prescrip- 
tion for the contour of In P/P' = Ht to be in the upper 
complex plane. In this we recover the fact that when the 
anti-commutator in the vacuum is P.V.(l/t), the corre- 
sponding vacuum Wightman function is l/(t — it). 

In Fig. 5, we plot the ratio T^/Tqh as a function of oj 
for various values of A, and for = 0. We first observe 
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FIG. 5. Plot of the ratio TL/Tgh as a function of cj/Tgh 
for various values of A. We work with a massless field 
with g — gcrit, for a detector localized in the center of 
the patch (X — 0), and with Ty, = 0. The values of A 
are 1 (continuous), 3 (dot-dashed), 5 (dashed), and 10 
(dotted.) 
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FIG. 6. Plot of To/Tgh — 1 in logarithmic scales as a 
function of A, where To is the low frequency temperature 
of massless fields with g = (? cr it and when t? = (ip- 
vacuum). We have also represented by a dotted curve 
the same quantity evaluated without dissipation when the 
state is the Bunch-Davies vacuum. 



that T u is constant for all frequencies from zero to a few 
multiples of Tqh- Hence, the Planckian character of the 
state is, to a high accuracy, preserved by dissipation, as 
was found in the presence of dispersion [15, 38, 39]. For 
higher frequencies, i.e., w/Tqh > 4, we were not able to 
study with sufficient accuracy because of the numer- 
ical noise associated to < 0.01. As in the dispersive 
case, we expect that the temperature function T w is mod- 
ified for uj > A. 

Secondly, when A is smaller than 5, i.e., when dissipa- 
tion is strong, we observe that the temperature is signifi- 
cantly (more than 5%) larger than Tqh- These deviations 
are further studied in Fig. 6, where we plot the deviations 
of To, the low frequency effective temperature, with re- 
spect to Tqh as a function of A. We observe that the 
deviation due to dissipation asymptotically follows 



^-1 ~ 



(6A)" 



(46) 



This law has been verified up to A = 10 3 . It has to be 
compared with the deviation due to quartic dispersion 
studied in Ref. [15]. This deviation is represented by 
the dotted curve, and scales as Tp isp /Xgh — 1 ~ e~ nX / 4 . 
In other words, the deviation due to (quadratic) dissi- 
pation decreases much slower than that due to (quartic) 
superluminal dispersion. The important lesson for black 
hole thermodynamical laws is that ultraviolet dispersion 
and dissipation both destroy the thermality of the state. 
This lends support to the claim that Lorentz invariancc 
is somehow necessary for these laws to be satisfied. 



relativistic out basis at fixed k to read out the state of 
(j). Alternatively, one can also use an out basis formed 
with stationary modes with fixed frequency lo. Indeed, 
at fixed ui, the momentum P u ~ |^/A| — > at large 
\X\, and dispersive effects are negligible. Hence ^(X), 
the stationary component of the field operator, decou- 
ples from the environment at large \X\, and can be an- 
alyzed using relativistic modes. As we shall see, this 
new out basis is not trivially related to the homogeneous 
one used in Sec. Ill because it encodes thermal effects 
at the Gibbons-Hawking temperature. Hence the covari- 
ance matrix of the new out operators will depend on rif. 
and Cfc of Eq. (35), but also on these thermal effects. At 
this point we need to explain why we are interested in 
expressing in a different basis a state which is fully char- 
acterized by rife and c^. The main reason comes from 
black hole physics. As shall be discussed in the next 
section, when certain conditions are met, the results of 
this section apply to the Hawking radiation emitted by 
dissipative fields. 

To compute the covariance matrix in the new basis, we 
recall some properties of the relativistic massless field in 
de Sitter. First, because of conformal invariance, the field 
operator splits into two sectors which do not mix, one 
for the right moving [/-modes with k > 0, and the other 
for the left moving V-modes with k < 0. In addition, in 
de Sitter, the time-dependence of all homogeneous modes 
can be expressed through <p(P) of Eq. (27), which here 
reduces to 



tp(P) = e lP/H /V2H , 



(47) 



B. Asymptotic correlations among right movers 

As explained in Section III A, at late time, the (j> field 
decouples from its environment. This allows to use the 



where P > 0. This mode has a unit positive Klein- 
Gordon norm, as can be verified using the Wronskian 
condition of Eq. (28). 

We introduce an intermediate basis constructed with 
the stationary "Unruh" modes ip u [40]. In the P- 
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representation, they can be written as [41] 

tp u = {P/H)- lu ' H - 1 x tp{P) . 



(48) 



They form an orthonormal and complete mode basis if 
u) g] — oo,oo[. The spatial behavior of the [/-modes is 
given by 



dP 



iPX 



HV2ir 



(49) 



We now introduce the alternative out basis formed of 
stationary modes which are localized on either side of 
the horizons, henceforth called R and L modes. They 
behave as Rindler modes in Minkowski space. For U- 
modes, the horizon is located at HX = —1, and these 
modes are 



Xu' R (X) = 0(1 + HX 



(xu u ' L (x)r = e(-i-Hx 



(1 + HX) iu / H 

(-1 - Hxy 



,/H 



(50) 



2uj 



where uj > 0. The first has a positive norm, while the 
second has a negative one. They are easily related to the 
Unruh mode by computing Eq. (49). Indeed, for uj > 0, 
one gets 



Vv = ^ Xu + Pu> (Xu ) 



(51) 



where coefficients a^j and ($~f are the standard 
Bogoliubov coefficients leading to the Gibbons-Hawking 
temperature H/2-k. They obey \PS l a S\ = e~™/ H . 
Asymptotically in the future and in space, the /7-part 
of the field operator can thus be expressed as 



0u (x) 



dk{a k e llcz <p k (t) + h.c.} 



d<jj{t 
d<jj{t 



rt(X) + h.c.} 



U,R ' 



(52a) 
(52b) 
(52c) 



W 1 



(Xu!(X)T + h.c} 



The U-part possesses a similar decomposition, and the V- 
modes are obtained from the [/-ones by replacing X —>■ 
—X, and R —> L. The xv modes are thus defined on 
either side of HX = 1. 

Using the above equations, the Unruh and the Rindler- 
like operators of frequency |w| are related by 



a U,L 
a V.L 
V.B, ' 










V o 










\ 








X 


a U 
ay 




PS* 


PS 


*S*) 







(53) 



We considered both U and V modes because our aim is 
to compute the covariance matrix of the R and L oper- 
ators in terms of n k and c k of Eq. (35), where c k mixes 



U and V modes. To do so, we first compute the covari- 
ance matrix of the Unruh operators. When working with 
states that are invariant under the affine group, n k and 
Cfc of Eq. (35) are independent of k. This implies that the 
covariance matrix of the Unruh operators is independent 
of uj. Indeed, using 



°° dk 



iui/H -1/2 



(54) 



which follows from the Fourier transforms Eqs. (52a) 
and (52b), one verifies that the independence of k 
implies that of w. As a result, introducing Uj — 

, a^ w , ay", dy^^j , the covariance matrix of Unruh 

operators reads 



C = Tr 



5(uj — uj') x 



/ n k Cfc 

n k c* k 

Cfc n k 

V 4 n k 



(55) 



where n k and c k are given in Eq. (23). 

Using the matrix of Eq. (53), and dropping the 
trivial factor of 5(u) — uj'), the covariance matrix of R 
and L operators is 



CS^B^CBl 
( 

= 2 



n u 
c* 



MV - 



uv 



rn 



MV* 



c* 



UV 



(56) 



where 



2n L 



+ 1=(K| 2 + |A?| 2 ) (2«fc 
c u = aS{PST^n k + l) , 
muj = 2Re(c k aSl3S) 



1) 



2c u J = (aSfc k + [((3Sfc k 



(57a) 

(57b) 
(57c) 
(57d) 



The first two coefficients concern separately either the U, 
or the U-modes. They fix the spectrum and the strength 
of the correlations. The last two concern the U — V mode 
mixing, and are proportional to c k . 

Considering the coherence amongst pairs of [/-quanta, 
i.e., ignoring the V- modes, as in Eq. (25), we define 



1 



|cj 2 /n u 



Using Eq. (57), we obtain 



A" — 
°u — n 



H\2 



n k (n k - 

+ \PS\ 2 ) 



-A 

n k 



\P 



H\2 



(58) 



(59) 



We see that 8u does not depend on c k . This is to be 
expected since c k characterizes the correlation between 
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FIG. 7. Figure for 8^ with uj = H for massless field with 
critical coupling g — g C rit- In the infalling vacuum, for — 0, 
the non-separability found for the massless relativistic case is 
preserved as long asA/ff = A>l/5. When the environment 
is characterized by a temperature ^ 0, the entanglement 
is preserved as long as T^, < V HA/2, as explained in the text. 

modes of opposite momenta, and since there is no U — V 
mode mixing for two dimensional massless fields. More 
importantly, Eq. (59) is valid irrespectively of the tem- 
perature of the environment T^. We can thus study how 
the separability of [/-quanta is affected by T^. The cri- 
terion of non-separability, 6y < 1, gives 

\B H \ 2 = - > n K T 1>) ( 60 ) 

where rik{T^) is plotted in Fig. 2. Using this Figure, in 
Fig. 7 we study In 8% with uj = H as a function of A 
and "d = T^/H. At zero temperature = 0, we see 
that the pair of [/-quanta with ui = H is non separable 
for A > 0.2, i.e., for a rather strong dissipation since 
A = H/5. Using Eq. (60) we see that this is also true for 
all quanta with uj/H < 1. More surprisingly, when A is 
high enough, this pair is non-separable even when > 
Tqh, i.e., when the environment possesses a temperature 
higher than the Gibbons-Hawking temperature. Indeed, 
whenever < v HA/2, the pair is non-separable, as all 
pairs with smaller frequency lo. 

In other words the quantum entanglement of the low 
frequency [/-pairs of quanta is extremely robust when 
working with dissipative fields which are relativistic in 
the infrared. The robustness essentially follows from the 
kinematical character of the transformation of Eq. (53) 
which relates two relativistic mode bases. It is also due to 
the fact that n k , the number of U—V pairs created by the 
cosmological expansion, remains negligible in Eq. (57) as 
long as 1 < A/H, and < T GH (A/H)^ 2 . 



V. BLACK HOLE RADIATION 

We now explain when and why the above results apply 
to the Hawking radiation emitted by dissipative fields. 
We shall be more qualitative than in the former sections 
because several approximations are involved in the corre- 
spondence between de Sitter and the black hole case. Our 
main aim is to establish that the spectrum of Hawking 
radiation, and the associated long distance correlations 
across the horizon, are both robust when dissipation oc- 
curs at sufficiently high energy with respect to the surface 
gravity, as was anticipated in Refs. [20, 22]. 

The robustness shall be established by studying the 
anti-commutator of Eq. (A3b), and showing that its 
asymptotic behavior is governed by Eqs. (57a) and (57b). 

Firstly, being covariant, the action of Eq. (2) applies 
as such to any black hole metric endowed with a pre- 
ferred frame described by a time-like field u. 2 Secondly, 
the correspondence with de Sitter becomes more precise 
when working with stationary settings. At the level of 
the background, this means that there is a Killing field 
K t , and that u commutes with K t . In this case, the 
metric can be written as 

ds 2 = -dt 2 + (dX - v(X)dt) 2 . (61) 

As in Eq. (18), t, X are defined by dt = vJ^dx^, and dx = 
s ffdtJ,i where u is a stationary and freely falling unit 
time-like field. In the present case, it is no longer unique 
because the system is no longer translation invariant. It 
belongs to a one parameter family, where the parameter 
can be taken to be the value of v at spatial infinity [44]. 
When the preferred field u is freely falling (as we shall 
assume for simplicity), this residual invariance is lifted 
by working with u s = u. 

By stationary settings, we also meant that the state 
of the environment is stationary. This implies that the 
noise kernel of Eq. (16) only depends on t — t' when eval- 
uated at X,X' , along the orbits of the Killing field K t . 
When these stationary conditions are met, the (driven 
part of the) anti-commutator of <f> is (exactly) given by 
Eq. (A3b), where the two kernels G% t and N 1 ^ are now 
defined in the black hole metric of Eq. (61). 

As a result, to compare the expressions of G% C (X, X') 
evaluated in de Sitter and in Eq. (61), it is sufficient to 
study Gf et and A". To establish the correspondence with 
controlled approximations, four conditions are necessary: 

• the state of i/j should be the same, 



2 When completing this work, we became aware of Refs. [42, 43] 
where similar issues are discussed. While the model is similar to 
that of Eq. (2), the preferred frame is taken at rest with respect 
to the orbits of the stationary Killing field K t . This means that u 
is space-like in the supersonic region. Unlike what is claimed, we 
believe that dissipation will necessarily engender an instability. 
More generally, we have not been able to follow the mathematical 
developments of these works. 
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• the black hole surface gravity n = H, 

• the near horizon region should be large enough, 

• the dispersive and dissipative scales should be both 
much larger than k. 

The first condition is rather obvious and needs no justifi- 
cation. The second and the third conditions concern the 
metric and the field u. To characterize the near horizon 
region (NHR) explicitly, we shall use 



v = -l + Dtanh(KX/D) 
1 + nX + D 0(kX/D) 3 



(62) 



which possesses a future (black hole) Killing horizon at 
X = 0. The NHR is defined by the region \kX\ < D/2 
where v is approximately linear. Hence it is a portion 
of de Sitter space with H = k, see Eq. (18). It should 
be emphasized that the mapping also applies to the u 
field. In fact, when u is freely falling, the only scalar 
quantity which is involved in the mapping is its expansion 
evaluated at the horizon: Qq = — V M u M = k. Hence, in 
the NHR, the orbits of u coincide with those found in 
de Sitter. (When u is accelerating, both Oo and the 
acceleration 70 must match, see Eq. (A10) and footnote 
4 in Ref. [15].) 

Using Eq. (62), the third condition means that D can- 
not be too small. This condition was found in Ref. [38] 
when considering the spectral deviations of Hawking ra- 
diation which are due to high frequency dispersion, see 
also Refs. [39, 45, 46]. For quartic dispersion, these devi- 
ations are small when D 3 ^ 2 ^> k/A. In this case, the non- 
trivial dispersive effects all occur deep inside the NHR, 
i.e., in a portion of de Sitter space. Moreover, at fixed 
k/A, the spectral deviations increase when D decreases. 
We shall see below that these facts also apply to dissipa- 
tive fields when the above four conditions are met. 



A. The stationary noise kernel 

When considering the model of Eq. (2) in the metric 
Eq. (61) with u freely falling, the noise kernel of 
Eq. (A3b) is 

N"(X U X 2 ) = 7 (-di)7(-&) (-UJ + V^diV^) (63) 
(iuj + ^282^2) JdqG^(X 1 ,X 2 ,q) , 

where Vi = v(Xi) and di = dxi- The stationary ker- 
nel of the last line is the Fourier transform of the anti- 
commutator of -0, see Eq. (36). To compute it we use the 
fact that the factor a(r, z) of Eq. (6) is now given by (see 
Eq. (55) in Ref. [22]) 



a(X,t) = v(X)/v{z{X,t)) 



(64) 



As in Eq. (6), z labels the orbits of u. It is here com- 
pletely fixed by the condition that z = X when t = 0. 



Since the orbits are solutions of dX/dt = v, z is implicitly 
given by 



x 



dXi 



t 



(65) 



Using the above equations to re-express the 5(z — z') of 
Eq. (36), one finds 



5{At-f£dX/v) 



G ajC ^(At,X 1 ,X 2 ;q) = 



2n a 



(66) 



COS (£l q At) 



Its Fourier component with respect to At is trivially 



G^ c< AXi,X 2 ;q) — 



juAt 12 2n q + l 



a, 



cos(fi ? Ai 12 ) , (67) 



where At 12 = f£ dX / v is the lapse of time from X2 to 
X\ following an orbit z — est. which connects these two 
points. Since the settings are stationary, these orbits are 
all the same, as can be seen in Fig. 9. 

Using Eq. (63), the noise kernel is explicitly given by 



N"(X U X 2 ) =7(-di) 7 (-d 2 ) 



jdq(2r 



ibjAt 12 



(68) 



1) f2 9 cos {p, q At\<x) 



This kernel is local in that it only depends on g M1/ and 
w M between X\ and X2- Hence, when evaluated in the 
black hole NHR, it agrees, as an identity, with the cor- 
responding expression evaluated in de Sitter. 

In conclusion, we notice that this identity follows from 
our choice of the action of Eq. (2) . Had we used a more 
complicated environment, this identity would have been 
replaced by an approximative correspondence. In that 
case, the correspondence would have still been accurate 
if the propagation of ip had been adiabatic. As usual, 
this condition is satisfied when the degrees of freedom of 
ip are "heavy", i.e., when their frequency Q q ~ A 3> K. 



B. The stationary Gf ct 

The stationary function G" et (X,Xi) obeys Eq. (12), 
which is a fourth order equation in dx when working with 
Eq. (39). Depending on the position of X and X\, its 
behavior should be analyzed using different techniques. 
Far away from the horizon, the propagation is well de- 
scribed by WKB techniques since the gradient of v is 
small. Close to the horizon instead, the WKB approxi- 
mation fails, as in dispersive theories [45]. In this region, 
the P-representation accurately describes the field prop- 
agation, and is essentially the same as that taking place 
in de Sitter. Therefore, the calculation of G% C (X,X') of 
Eq. (A3b) at large distances boils down to connecting the 
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de Sitter-like outcome at high P to the low momentum 
WKB modes. As in the case of dispersive fields, the con- 
nection entails an inverse Fourier transform from P to 
A-space in the intermediate region 77, see Fig. 9, where 
both descriptions are valid [20, 45, 47-49]. In the present 
case, these steps are performed at the level of the 2-point 
function rather than being applied to stationary modes. 
In fact, we shall compute G ac through 

G^(X,X') = JJ™ dP 1 dP 2 G? et (X,P 1 ) 

xG? et *(A',P 2 )A"(Pi,P 2 ) , 

where the two G" ct are expressed in a mixed X, P- 
representation. The early configurations in interaction 
with the environment are described in P-space, while the 
large distance behavior is expressed in A-space. 

Let us give here only the essential points, more details 
are given in Appendix C. The validity of the whole pro- 
cedure relies on a combination of the third and the fourth 
condition given above, namely max(l,P~ 2 ) <C A/k, and 
is limited to moderate frequencies, i.e., < U) ~ k <C A. 

For simplicity, we consider massless fields. Then 
A/it » 1 guarantees that the infalling V^-modes essen- 
tially decouple from the outgoing [/-modes because the 
only source of U — V mixing comes from the ultraviolet 
sector. Hence, at leading order in k/A, it is legitimate 
to consider only the [/-modes. For massive fields with 
m <C A, the discussion is more elaborate but the main 
conclusion is the same: the properties of the Hawking 
radiation are robust. 

For massless fields, at fixed uj, the propagation of the 
[/-modes is governed by the effective dispersion relation, 
see Eq. (34), 

= u> - v(X)P = \/F 2 - T 2 . (70) 

As long as P <C A, the [/-sector of G^ et behaves as for a 
relativistic field, since %/P 2 - T 2 - P(l + 0(P/A)). In- 
stead, when P > A, the dispersive and dissipative terms 
weighted by / and T cannot be neglected in Eq. (12). To 
characterize the transition from these two regimes, we 
consider the optical depth of Eq. (32). When working at 
fixed w, one finds 

WXl) = /p ^ P> dxv[x\p>)] ' 

where X U (P) is the root of Eq. (70), as is P bJ {X) when 
using X as the variable. The first expression governs 
G rot in the NHR where dxv ~ K is almost constant, see 
Eq. (31). To leading order in T/P <C 1, which is satisfied 
everywhere but very close to the horizon, the second ex- 
pression governs G rc t in A-space. Since Vg V — l/d u P is 

the group velocity in the rest frame, 1^ = J t dt'T(P u ), 



where the integral is evaluated along the classical outgo- 
ing trajectory. It should be noticed that, when consid- 
ered in A-space, X u applies on the right and the left of 
the horizon. In the R region, v gT > 0, while it is negative 
in L, so that in both cases X w > when Pi > P > 0, i.e., 
when Pi is in the past of P. 

To characterize the retarded Green functions of 
Eq. (69), we compute X u in the mixed representation, 
in the limit where Pi is large enough so that X ul (Pi) 
is deep inside the NHR, while X is far away from that 
region. For simplicity, we consider the case of Eq. (39) 
with g = <7 cr it- In this case, only the dissipative effects 
are significant, 3 and one finds 

W-lk + ^TF' <72 > 

where vr (vl) is the asymptotic velocity on the right 
(left) side. From the second term, we learn that \kX\ 
should be much smaller than A/k for the Hawking quanta 
not to be dissipated. Since we work in the regime A/k> 
1, this condition is easily satisfied. We notice that a 
similar type of weak damping effect of outgoing modes 
has been observed in experiments [50] . 

From the first term, we learn that 1^ gives an upper 
bound to the domain of P which significantly contributes 
to Eq. (69), namely P 2 < Art, as in de Sitter. A lower 
bound of this domain is provided by the 7 factors of 
Eq. (38). Using this equation and Eq. (A8), the integrand 
of Eq. (69) scales as 

T(P) ocPr(P)e- 2I ^ p ) cx P 3 e ^ 2 /AK . (73) 

Its behavior represented in Fig. 8. Hence, the relevant 
domain of P, i.e., when T is larger than 10% of its max- 
imum value, scales as 

0.36^ = P min < P < P max = 2.47kA . (74) 

Considered in space-time, since P ~ e _Kt , this limits 
the lapse of time during which the coupling to tp occurs. 
Interestingly, this lapse is given by nAt « 2, i.e., two e- 
folds, irrespectively of the value of A/k, and that of lu. It 
should be also stressed that nothing precise can be said 
about the domain of X which significantly contributes 
because the A- WKB fails when P is so large. One can 
simply say that it is roughly characterized by the interval 
[-A trans , A trans ], where A trans = A w=K (P min ) is given by 

rtXtrans ~ 3a/k/A . (75) 

This value defines the central region III, see Fig. 8 and 



3 In the case where j 2 « 1, dispersive effects are important and 
may limit the role of dissipation in the NHR. In that case, the 
decaying part of the field will contribute to Eq. (14). This sit- 
uation corresponds to what is found in the surface wave experi- 
ments [50, 51]. 
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FIG. 8. As a function of log 10 (P/«:), in a dotted line we 
plot exp{— X U1 (X, P)}, the optical depth of Eq. (71), evaluated 
for oj/k = 1, A/k = 400, and nX = 20. The solid curve 
represents T(P) of Eq. (73) for the same values, and D = 
0.99. The left dash-dotted line corresponds to the limit of the 
NHR: kX = D/2, here re-expressed as P = 2n/D. For lower 
P, in region 7, the de Sitter-like P-representation fails. The 
right vertical line indicates the upper limit of the X-WKB 
approximation, see Appendix C. For the adopted values, the 
region 77 where the P and the X descriptions are both valid 
has a finite size. We also see that T vanishes in region 7. 

Fig. 9. Using the profile of Eq. (62), A" trans is situated 
deep inside the NHR when k/uj^^ <C 1, where the criti- 
cal frequency w^ax i s given by 

<ax = A/J 2 . (76) 

Hence, when k/ui^H. <C 1, the coupling between 4> and 
if) is accurately described in the P-representation, and 
takes place in a portion of de Sitter. In addition, the 
connection between the high and low momentum propa- 
gation can be safely done in the intermediate region II, 
defined by K|A" trans | <C k\X\ < D, see Fig. 9, where, on 
the one hand, one is still in a de Sitter-like space since 
v is still linear in X, and, on the other hand, the low 
momentum modes can be already well approximated by 
their WKB expressions. Notice finally that this reason- 
ing only applies for frequencies oj <C wjjj. Indeed, when 
w = wjjjax, dissipation occurs around nX ~ D, i.e. no 
longer in a de Sitter like background. 

These steps are sufficient to establish that the results 
of Section IV B apply for cj -C w*^. In particular, 
Eq. (57a) implies that the spectrum of radiation is ro- 
bust (when the temperature of the environment is low 
enough, see Fig. 7). Namely, to leading order in k/A, 
the mean occupation number n u of quanta received far 
away is given by the Planck distribution at the standard 
relativistic temperature Th = k/2tt. As in dispersive 
settings, the real difficulty is to evaluate the spectral de- 
viations. In this respect, we conjecture that the leading 
deviations due to dissipation will be suppressed by pow- 
ers of /V w max- That is, they will be governed by the 
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FIG. 9. Null outgoing geodesies (dashed lines) on either side 
of the horizon at X = 0, and freely falling orbits z — est. (dot- 
ted) in the t, X coordinates of Eq. (61). As explained in the 
text, the nearby geodesies schematically indicate the space- 
time region where G ac (t, X, ti, Xi) is non- vanishing, when 
Kti — 2.5, and kX\ = 1.5, see Fig. 1 of Ref. [41] for the 
relativistic case. The two thick solid lines represent the re- 
gion where the noise kernel contributes to G ac , see Eq. (A3b) 
and Eq. (69). In the central region 777, the propagation is well 
described in P-space, and resembles to that found in de Sitter. 



composite ultraviolet scale of Eq. (76) which depends on 
the high energy physics, here with T quadratic in P, and 
on the extension D of the black hole NHR. This second 
dependence is highly relevant when Del. 

Together with the robustness of the spectrum, one 
also has that of the long distance correlations across 
the horizon between the Hawking quanta and their part- 
ners. These correlations are fixed by the coefficient c u of 
Eq. (57b). To get the space-time properties of the pat- 
tern, one should integrate over u>, i.e., perform the inverse 
Fourier transform of Eq. (19b), because it is this integral 
that introduces the space-time coherence [20, 41, 52]. 
In Fig. 9, we have schematically represented the anti- 
commutator G ac (t— ti, X, X{) in the t—ti, X plane when 
X\ is taken far away from the horizon. 

• Far away from the NHR, in regions In and II, for 
D < \kX\, the characteristics of the field follow null 
geodesies, see Eq. (C10). Since v ~ est., they no 
longer separate from each other. Hence, at large 
distances, the space-time pattern obtained by fix- 
ing one point [52], and the equal time correlation 
pattern [53] , will be the same as those predicted by 
a relativistic treatment. 
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• In the two intermediate regions II r and II l, for 
-^trans < \ K X\ ?5 D, the characteristics separate 
from each other following 5X ~ e Kt since their be- 
havior is already close to the relativistic one. This 
pattern is obtained by considering two-point func- 
tions with one point fixed, or wave-packets [20]. It 
is interesting to notice that it cannot be obtained 
by considering equal time correlations, since these 
develop only outside the NHR, for \kX\ > D [41]. 
Indeed as long as X and X' are in the NHR, the 
(approximate) de Sitter invariance under K z , see 
Appendix A, implies that J dujG^ c only depends 
on X - X 1 . 4 

• In the central region III is the region where the 
configurations of the </> field are driven by the noise 
kernel. In Fig. 9 the two thick solid lines indicate 
the space-time locus where the interactions involv- 
ing the configurations selected by t\ , X\ are taking 
place. 5 In this central region 777, the propagation is 
well described in P-space, and corresponds to that 
found in de Sitter, see Eq. (C14) and Eq. (C15). 

In brief, when K/^raX ^ 1 an d w / w mS ^ 1) tne 
non-trivial propagation only occurs deep inside the NHR 
which is a portion of de Sitter space. This implies that 
n u and c w are, to a good approximation, given by their 
de Sitter expressions of Eq. (57). Given that these (exact) 
expressions hardly differ from the relativistic ones when 
k/A <C 1, we can predict that, when computed in a black 
hole metric, these two observables are robust whenever 
the finiteness of the NHR introduces small deviations 
with respect to the de Sitter case. For cj/lu^ x <C 1, 
this is guaranteed by ^/^mS <C 1. 



VI. CONCLUSIONS 

In this paper we used the dissipative model of Ref. [22] 
to compute the spectral properties and the correlations 
of pairs produced in an expanding de Sitter space. The 
terms encoding dissipation in Eq. (2) break the (local) 
Lorentz invariance in the ultraviolet sector. Yet, they are 



4 From this observation, we learn that the correspondence between 
the physics in black hole metrics and in de Sitter is not merely a 
convenient way to obtain riu and c u in Appendix C 2. It actually 
shows up in the NHR when computing observables, such as the 
mean value or the 2-point correlation of p = u^u^T^y. Moreover, 
it ceases when leaving this region. In this sense, the Hawking 
effect only develops, or separates, from its de Sitterian roots for 
\kX\ > D, and furthermore, this separation is adiabatic. 

5 To get Fig. 9, we have filtered out low frequencies. This amounts 
to consider a wave packet rather than the two-point function, see 
Eqs. (38,39) in Ref. [41]. Had we considered G ac , the thick lines 
would have extended back in time, because the coupling between 
ip and 0^ is centered around nt^ = lno; + est., which fixes the 
blue-shift for P u a w to reach V kA, see Eq. (74). 



introduced in a covariant manner by using a unit time- 
like vector field u which specifies the preferred frame. In 
addition, the unitarity of the theory is preserved by cou- 
pling the radiation field (f> to an environmental field ip 
composed of a dense set of degrees of freedom taken, for 
simplicity, at rest with respect to the u field. Again for 
simplicity, the action is quadratic in <^>, , 0, and the spec- 
tral density of -0 modes is such that the (exact) retarded 
Green function of (f> obeys a local differential equation, 
see Eq. (10) and Eq. (12). 

By exploiting the homogeneous character of the set- 
tings, we expressed the final occupation number rife, and 
the pair-correlation amplitude Cfe, in terms of the noise 
kernel TV and the retarded Green function, see Eq. (35). 
Rather than working with integrals over time as usu- 
ally done, we used the proper momentum P = k/a(t) to 
parametrize the evolution of field configurations. Hence, 
Eq. (35) can be viewed as flow equations in physical mo- 
mentum space. This possibility is specific to the residual 
symmetry group found in de Sitter space when the u field 
commutes with the two Killing fields K t and K z . These 
group theoretical aspects are explained in Appendix A. 
The key equations are Eq. (A6) and Eq. (A8) which show 
how the P-representation is related to the invariant dis- 
tances, to the homogeneous representation of Eq. (20), 
and to the stationary one. This representation is ex- 
tended to Feynman rules and Schwinger-Dyson equations 
of (relativistic) interacting field theories in Ref. [54]. 

We numerically computed and Cfe in Section III. 
When considering a massless field, rife and the strength 
of the correlations are plotted as functions of the scale 
separation A/H, and the temperature of the environment 
T^/H, in Fig. 2. The robustness of the relativistic results 
is established in the limit of a large ratio A/H. The key 
result concerns the threshold values of the parameters, 
see the locus Sk = 1 on the right panel, for which the 
final state remains non-separable, i.e., so entangled that 
it cannot be described by a stochastic ensemble. Vari- 
ous criteria of non-classicality are compared in Appendix 
B. This analysis was then extended to massive fields, see 
Fig. 3, and to the consequences of varying the relative im- 
portance of dissipative and dispersive effects, see Fig. 4. 
As expected, the quantum coherence is lost at high cou- 
pling, and when the temperature of the environment is 
high enough. 

In Section IV we exploited the stationarity, and we 
studied how the thermal distribution characterizing the 
Gibbons-Hawking effect is affected by dissipation. As in 
the case of dispersion [15], we found that the thermal 
character is, to leading order, robust. We also computed 
the deviations of the effective temperature with respect 
to the standard one Tqh = H/^tt, see Fig. 5 and Fig. 6. 
In preparation for the analysis of the Hawking effect, we 
studied the strength of the asymptotic correlations across 
the Killing horizon between (right) moving quanta with 
opposite frequency. Quite remarkably, we found that the 
pairs remain entangled (the two-mode state remains non- 
separable) even for an environment temperature exceed- 
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ing T G h = H/2-k, see Fig 7. 

Finally, in Section V we extended our analysis to black 
hole metrics. When four conditions are met, we showed 
that the above analysis performed in de Sitter applies to 
Hawking radiation. The inequality which ensures the va- 
lidity of this correspondence is k/oj^H- <C 1, where w^ 1 ^ 
is the composite ultraviolet scale of Eq. (76). It depends 
on both the microscopic scale A, and D, which fixes the 
extension of the black hole near horizon region where the 
metric and the field u can be mapped into de Sitter. The 
validity of the correspondence in turn guarantees that, to 
leading order, the Hawking predictions are robust - even 
if the early propagation completely differs from the rela- 
tivistic one, see Fig. 9. This establishes that when leaving 
the very high momentum P ~ A (trans-Planckian) region 
and starting to propagate freely, the outgoing configura- 
tions are "born" in their Unruh vacuum state [34, 55, 56]. 
The microscopic implementation of this state in dissipa- 
tive theories is shown in Eq. (C16). 

To conclude, we point out that the correspondence be- 
tween black hole physics and de Sitter physics tells us 
that, as in the case of dispersive theories [38, 45], the 
leading deviations with respect to the relativistic expres- 
sions should be suppressed as powers of K/^maxi i-e., they 
should be governed by the extension of the black hole 
NHR which is a portion of de Sitter space. 
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Appendix A: AfHne group and P-representation 

We remind the reader that the affine group is the sub- 
group of the de Sitter isometry group which is generated 
by the Killing fields K z = d z \t and K t = dt\x, which pos- 
sess the following commutator [K z) K t ] — —HK Z . The 
definition of the coordinates t,z,X is given in Eq. (17) 
and Eq. (18). In de Sitter space, there are two geomet- 
rical invariants under this group. Using the coordinates 
t, X, they read 

Ai = e ff (*-*') , (Ala) 
A 2 = le^^ 1 ''/ 2 - X'e H ^^l 2 . (Alb) 



They are linked to the de Sitter invariant distance by 

A 2 = A 2 (A a - -^) 2 . (A2) 

The distances Ai, A 2 can be also defined in a coordinate 
invariant manner. The interested reader will find the 
expressions at the end of this Appendix. 

When working with states that are invariant under the 
affine group, the n-point correlation functions only de- 
pend on Ai and A 2 evaluated between the various pairs 
of points. Hence, any 2-point functions G any (x, x') can be 
written as G a ny (Ai(x,x'), A 2 (x,x')). However, it turns 
out that it is not convenient to use Ai, A 2 to compute 
Eq. (15), and this even though the four integrals of that 
equation can be easily expressed in terms of two over Ai 
and two over A 2 . The reason is that the integrals over 
the A 2 are convolutions. Hence, it is appropriate to work 
with the Fourier transform with respect to A 2 because, in 
this representation, Eq. (15) contains only two integrals. 

The fact that only two variables are needed is not a 
surprise, given the homogeneity (stationarity) of the set- 
ting. Indeed using G k (t,t') (G U (X,X')) of Eq. (19), one 
immediately has 

Gt(t,t')=J JdhdhGUt^G^it'^N^t^h) , 

(A3a) 

G^(X,X')=f JdX 1 dX 2 G^ ct (X,X 1 ) (A3b) 

xG? e l(X',X 2 )N u (X x ,X 2 ) . 

To understand the relationship between these two re- 
presentations, it turns out that the most convenient vari- 
ables are the proper momenta P = and P' = 
|SffPu|- The reasons for this are many. Firstly, P is in- 
variantly defined; secondly, A x is easily expressed in P, P' 
space; thirdly, so is the variable conjugated to A 2 ; and 
fourthly, P can be attributed to the field itself, so that 
one can easily take the even (anti-commutator) and the 
odd part of the 2-point functions. Let us explain these 
reasons. 

Once the de Sitter group is broken in a way which 
preserves the affine group, P is invariantly defined as 
the momentum associated with the orthogonal fields 
us , sg which commute with K t and K z , and where us 
is geodesic. In our case, we work with the preferred field 
u = tig, but this needs not be the case for P to be un- 
ambiguously defined as P 2 — (sgP M ) 2 . 

Since P — ke~ Ht , Ai is simply 

A 1 (x,x')=P'/P>0. (A4) 

In addition, the momentum conjugated to A 2 , defined by 
P = (?a,|Au i s given by the geometrical mean 

P = P v / Ai= sgn(P)/PP 7 . (A5) 

The first equality follows from A 2 y/A^ = X + f(t, t' , X'), 
and P = dx\t,t',x'- The second one follows from 
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Eq. (A4). Hence, the Fourier transform of G any (Ai, A2) 
with respect to A2, 



Gany (P j P') 



(PP 



H 



dA2e -iVPP^ SS n( P )A 2 G anv ( -,A 2 



'P' 

'" y 1 p 



(A6) 



only depends on P and P'. Moreover, if one imposes the 
isotropy of the setting, G any (Ai, A2) is even in A2, and 
G any (— P, — P') = G any (P,P'). Hence, in this case, all 
the information is contained in G my (P, P'). 

The important point is that G any (P, P') defined by 
Eq. (A6) coincides with the l.h.s. of Eq. (20). In addition, 
starting with the stationary representation of Eq. (19b), 
one can also verify that the double Fourier transform 



G a V(P ! P0=/™e-' PV ^ - G; uv ,X. X' ) 



2n 



has automatically the following structure 



C any (P, P ) 



(P/P') 



-iu/H 



PP 1 



-G any (P, P') 



(A7) 



(A8) 



where G any (P,P') is given by Eq. (A6). Together with 
Eq. (20), Eq. (A6) and Eq. (A8) are the key equations of 
this appendix: Whenever a 2-point function G any (x, x') 
is invariant under the affine group, its Fourier transforms 
Gl uy {t,t') and G£ ny (P,P') are related to G any (P,P') of 
Eq. (A6) by Eq. (20) and Eq. (A8) respectively. Finally, 
the antisymmetry of the commutator G c is expressed as 
G C (P',P) = -G C (P,P')* while the symmetry of G ac 
gives G ac (P',P) = G a , c (P,P')*. 

To conclude this Appendix, we express Ai and A2 in 
covariant terms. The log of Ai is given by the line inte- 
gral of it ff from x to x', that is 



InAi 



-H 



(A9) 



This is an invariant expression. Indeed, on the one hand, 
since u s is geodesic, is an exact 1-form and the 

above integral does not depend on the path. On the other 
hand, u s is the only (time-like) unit geodesic field that 
commutes with K z and K t . Since, Eq. (A2) gives A2 as 
a combination of Ai and A which are both invariantly 
defined, so is A 2 . 6 



6 If one wishes, A2 can also be seen as the integral of s^dx^ , the 
1-form associated to the vector orthogonal to « ff . Since this form 
is not exact, one has to specify the contour from x to x'. Using 
the t, z coordinates, one should go at fixed z from t to (t + t')/2, 
then vary z at fixed time until ((f.+t')/2, z'), and vary t at fixed z 
until (t',z'). Any different contour would give some combination 
of Ai and A2. 



We notice that the preferred frame fields u, s have not 
been used. But, if one wishes, they can be used. Indeed 
any couple of orthogonal fields u, s which commute with 
K z and K t are related to u s , s s by 

u s = (9m + js) /H, s K = (6s - -yu) /H , (A10) 

where the constant expansion is = — V^u' 1 , and where 
the constant acceleration is 7" = u^V — 7s". 



Appendix B: Non-separability and Cauchy-Schwarz 
inequalities 

In this appendix, we consider homogeneous states. 
This implies that the state factorizes as 



P = 



IP2 fc) 



(Bl) 



fe>0 



(k) 

where p 2 fixes the state of the two-mode system k, — k. 
This also implies that rik and c& of Eq. (23) only depend 
on k. To be general, we work with ^ which 
means that the state is anisotropic. Our aim is to com- 
pare four inequalities relating the norm of Ck to rik and 
n_k which allow to distinguish quantum from classical 
correlations, for a recent review, see e.g., Ref [57] 



1. CS inequality in quantum mechanics 

Any quantum state (density matrix) p defines a (posi- 
tive) scalar product on operators by: 



(A, B) p = Tr (pA^B 



(B2) 

The corresponding Cauchy-Schwarz (CS) inequality im- 
plies 

Tr [pA^B^j 2 < Tr (pA^X) x Tr (p&B) , (B3) 

:n ap 
one gets 



When applied, to Eq. (Bl) with A — and B = a^L k , 



\c k \ < «k(n-k + l) ■ (B4) 
When «k = ft-kj on e obtains Eq. (24). 

2. Separability 

A bi-partite state is said separable [6, 58] when it can 
be written as 



_ V d o [k) 



(B5) 



where p n > 0, and where the two-mode states p 2 n are 

factorized p^^ = ® /5l" k ^ . The operators pi ±k ' ) are 
density matrices for each one-mode system at fixed k. 
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The structure of these states defines a new scalar prod- 
uct. It is given by 



{X, Y) scp = Pn Tr (p™x) Tr (p£> Y 



(B6) 



where X, Y are arbitrary operators. Considering oper- 
ators that act on one sector only, i.e., A — A <Ei 1 and 
B = 1 ® B, one finds 



Tr (ftepAtfi) = ^ 



(k) • -(-k) 



Tr (fteplti) = £„£„ Tr ($>A*A 



3i k) 



(A, A) 



{A, B) sep , 

(B7a) 

(B7b) 



scp ) 



where the quantities with a bar are the expectation values 
involving only one-mode states 



Cl k) = Tr 



(B8) 



The inequality in Eq. (B7b) comes from the positivity of 
Tr(pS^itin) applied to |„ = A - A^ , which gives 

2 



Tr 
Tr 



(p^a^a) 



)>- 


-r(k) 
A n 


)>- 


-j(k) 



(B9) 



The crucial point here is that the bound is insensitive 
to the ordering of A and A 1 * . Therefore, when applying 
the CS inequality associated with the scalar product of 
Eq. (B6), i.e., |(X,Y) sep | 2 < (X,X) scp x (Y,Y) 8ep , to 



X = fik and Y = a_ k , the strongest bound is 



\c k \ < n k n_ k 



(BIO) 



The only difference with Eq. (B4) is that n_k + l has been 
replaced by n_k by virtue of Eq. (B9). In conclusion, the 
inequalities of Eq. (26) characterize the quantum states 
which are non-separable. 



3. CS inequality for commuting variables 

When treating the statistical properties of radiation 
fields using classical (i.e., commuting) random variables, 
destruction and creation operators are replaced by com- 
plex variables ak — > ak, a k ~* a k' an( ^ the two-mode state 
is described by a two-mode distribution P^^ak, a_k), 
which is positive and normalized. The strength of the 
correlations is described by <?k,-k given by the following 
quartic expression: 



9u,-k = II d 2 a k d 2 a_ k P 2 (A:) (ak,a-k) |a k | 2 |a_ k 



(Bll) 



The spectra can be described by gk,k and <7_k,-k which 
are given by the above integral when the integrand is 
respectively |ak| 4 and |a_k| 4 - Eq. (Bll) can be seen as 
the scalar product on square integrable functions with 

(k) 

measure P 2 . The corresponding CS inequality is [59, 



|ffk,-k| < 5kk x ff-k-k 



(B12) 



When taking into account the quantum character of 
the radiation, 5k, -k of Eq. (Bll) should be replaced by 
the normal ordered expression [59] 

^ k = Tt(/£ fc) 4fiLM-k) • (B13) 

When considering Eq. (B3) with A = fik = a^ak and 
B = n_k, in the place of Eq. (B12), one gets 

2 



QM 
£/k,-k 



<«k+»k)x(ff- 



k.-k 



+ n_ k ) 



(B14) 



Hence, the exist quantum states for which M k violates 
Eq. (B12). 

We now relate Eq. (B12) to Eq. (BIO). We first no- 
tice that, in general, these inequalities are not related. 
However, Gaussianity implies that these statements are 
in fact equivalent. Using Wick's theorem, in Gaussian 
homogeneous states we get 



QM 
i/k,k' 



2n k 



™k«--k + \Ck\ 

Hence, Eq. (B12) reads 



if k = k' , 
if k = -k' 



n k n-u + \c k \ < 2n k n_ k 



(B15) 



(B16) 



which is equivalent to Eq. (BIO). Hence, when non- 
Gaussianities are known to be small, the violation of 
Eq. (B12) can be used as a criterion for the non- 
separability of the state. 

It is easy to find examples of separable (non-Gaussian) 
states which violate Eq. (B12). For instance, the two- 
particle factorized state 



„(fe) 
P\ = 



gives 



(al |0 k ) (0 k | a k ) <8 (al k |0_ k ) (0_ k | a. 



Q QM _ Q QM _ -, 
fk.k — u ' i/k,-k — 1 ' 



(B17) 
(B18) 



so that Eq. (B12) is violated. We are not aware of non- 
separable states which respect Eq. (B12). 



4. Subfiuctuant mode 

We show that non-separable states possess a subfiuc- 
tuant mode whose variance is smaller than that of the 
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vacuum. In the isotropic case, the proof can be found in 
Ref. [9]. Below, we extend the proof to the anisotropic 
case n k ^ ?i-k- 

To obtain the subfiuctuant mode, we diagonalize the 
2x2 covariance matrix Tt^^ {W" 1 , W}) with W — 
(a_k,ft k ) by a rotation, and not by a U(l, 1) transfor- 



mation (a Bogoliubov transformation). The operators 



L k = cos£e a_k + sin£e a k , 
Sk = - sin£e~ ie a_ k + cos£e ie a k , 



(B19) 



define the super- and the subfiuctuant mode, and the two 
angles are 



cos(2£) = (n k - n_ k )/y (n k - ?i- k ) 2 + \c k \ , (B20a) 
(9 = arg(c k )/2 . (B20b) 

One verifies that Tr(/5{5k, L^}) = 0, and that the 
spread of the subfiuctuant mode is 



Tr(/3{£k, = n k + n_ k + l- \j (n k - n_ k ) 2 + 4 \c k \ 2 . 

(B21) 

Using Eq. (BIO), one establishes that Tr(p{S" k , < I 
implies that the state is non-separable. QED. 



Appendix C: Flux and long distance correlations 

The expressions for the asymptotic flux and the corre- 
lation pattern are both encoded in Eq. (A3b). To obtain 
them, we need two things. Firstly, we need to character- 
ize G^ et from the asymptotic region down to the NHR. 
To this end, we should perform a WKB analysis of the 
stationary damped modes. Secondly, we need to connect 
the WKB modes with the high momentum de Sitter-like 
physics taking place very close to the horizon. 



1. WKB analysis 

At fixed w, using Eq. (10), Ddiss <^dec = implies that 
the decaying mode </)^ ec obeys 



(iw - 8 x v) (iu - vd x ) + F 2 (~d 2 x ) 
y(-d x )(iw - \/vdx\fv)i{dx) 



A 



(CI) 



,J dcc 



= 



The mode $j ec decays when displacing X along the di- 
rection of the group velocity. Hence, on the right of the 
horizon, the outgoing {/-mode decays when X increases, 
while it decreases for decreasing X < in the left region, 
see Fig 9. Hence, {/-modes spatially decay on both sides 
when leaving the horizon. 



As in the case of dispersive fields, we look for solutions 
of Eq. (CI) of the form 



<p(X) 



JJ X dX'Q„(X') 



(C2) 



where Qu>(X) is expanded in powers of the gradient of 
v(X). To first order, Eq. (CI) gives 



(u-v(X)Q u + iTf - (F 2 - T 2 ) = 

- \d x d Q [(a, - v(X)Q UJ + iT) 2 - (F 2 - T 2 



(C3) 



where the functions T > of Eq. (11) and F are evaluated 
for P = Q u . The leading order solution, the complex 
momentum Q^\x) = P^(X), contains no gradient, and 
obeys the complex Hamilton-Jacobi equation 



u-v(X)P + iT(P) = y f F 2 {P)-T 2 {P) = F(P) . (C4) 

As expected, this equation gives Eq. (1) since = oj—vP. 
To first order in the gradient, we get a total derivative 



Q 



i) 



-d x log 



(C5) 



Combining Eq. (C4) and Eq. (C5), we obtain the decay- 
ing WKB-mode 



VdecPO 



x e i!x dX'P.{X') 



IvgF(PS) 



(C6) 



To get this expression, we introduced Vg T = l/d^P^ 
which can be conceived as a complex group velocity. We 
also decomposed Pj 7 into its real part P u , and its imagi- 
nary part P(5- The oscillating exponential is the standard 
expression, while the decaying one is J dXP^. The lat- 
ter is equal to X u of Eq. (71) when working to first order 
in r/P, which is here a legitimate approximation. A 
preliminary analysis, similar to Eq. (A12) of Ref. [45], 
indicates that the corrections to Eq. (C6) are bounded 
A(i+v) 2 )■ Hence Eq. (C6) gives an accu- 
rate description everywhere but in the central region 777 
defined by KA trans of Eq. (75). 

Using Eq. (C6), the {/-mode contribution to the com- 
mutator is, for uj > 0, 

G-(X,X')=6(l UJ (X,X'))^ dcc (X) « rw (A'))* 

+ 9(l u (X',X))<p» grw (X) (<p% ec (X>))* , 

(C7) 

where the growing mode (fig rw satisfies Eq. (CI) with the 
opposite sign for the last term which encodes dissipation. 
The expression for u < is given by G^" = — (G^)* 
which follows from the imaginary character of G c in t, X 
space. We used the sign of I u in Eq. (C7) so that a 
similar expression is valid on the left of the horizon. 
Note also that Eq. (C7) cannot be used to estimate G" 
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across the horizon because the WKB approximation fails 
in region III. Note finally that Eq. (C7) is valid only for 
A\X-X'\ »1. 

Having characterized in quantitative terms the impact 
of dissipation, we now work in conditions such that the 
mode damping is negligible far away from this central 
region. That is, we work with X, X' obeying 



de Sitter for Pi > \J kA, when 



.diss 



of Eq. (76) obeys 



Xy 



« \x\ « yXT^ 



(C8) 



where the upper limit comes from the neglect of the sec- 
ond term in Eq. (72). Under these conditions, the anti- 
commutator of Eq. (A3b) is, for u) > 0, given by 

g: c (x,x>) =(2n w + 1) y&x) {v" R {x')y 

+ (<pl"(X))* ip - L "(X')] (C9) 
+ 2Re[ Ct ^(A)^(A')] , 

where and c w are constant because we are far from 
region III, and where the R and L out modes live on one 
side of the horizon and have unit norm. Being undamped, 
they are either relativistic, or, more generally, dispersive 
WKB modes. In the former case, they thus behave in 
the regions of interest, namely Ir/l and IIr/l, as 

viuj/k iuX/(l+v R ) 

^ R7 e{x)^^r.e{xy ^ 



2uj 1 



y/2u/(l + V R ) ' 



V L ' " V2u ' ^/\2lo/(1 + v l )\ 

(CIO) 

where Vril) is the asymptotic velocity in the region R 
(L, where 1 + vl < 0). As in de Sitter, the (positive unit 
norm) mode <p~y living in the L region has a negative 
Killing frequency. 

In Eq. (C9), n w and c w are unambiguously defined be- 
cause the R/L modes are normalized in regions Ir/l- 
Thus, they respectively define the spectrum emitted by 
the black hole, and the w-contribution of the correlation 
across the horizon. To compute them, we should find the 
equivalent of Eq. (35). To this end, we shall use Eq. (69), 
and exploit the fact that their values are fixed in the do- 
main of P given in Eq. (74). 



2. Connection with de Sitter physics 



ft/u^S <C 1. This second condition means that the high 
Pi behavior is governed by Eq. (31) and Eq. (A8). 

For simplicity we consider the massless case of Eq. (39) , 
when g 2 = 2. In this model, in de Sitter, using the Unruh 
modes of Eq. (48), the [/-mode contribution is 



G^ dS (X,P) = e- x ° 



^(X) (fl(P)^(P))* 

-(^(X))*(0(-P)^(P)) , 

(C12) 



where Iq is given in Eq. (32), and where we replaced its 
lower value P LU (X) <C VkA by because X is taken suffi- 
ciently large. Using Eq. (51), we can re-express Eq. (C12) 
in the R/L out mode basis. For w > we get 



G- dS (X,P) = e- I o 



x£W(x£(P)r 



(C13) 



In this we recover that the commutator possesses the 
same expression if one uses the in (Unruh) or the out 
mode basis. 

Eq. (C13) applies as such to the black hole metric in the 
regions IIr/l, kX 

trans -D/2, because G c bh 

obeys the same equations, and its normalization is fixed 
by the equal time commutators. In fact, in these re- 
gions the normalized black hole modes <p R , Lp~y coincide 
with the modes XwXT 1 of Eq- (50). Then, the WKB 
character of <p R , tpjy guarantees that Eq. (C13) applies 
further away from the horizon, in the regions defined by 
Eq. (C8). Hence, in these regions, we have 



G- BH (X,P) = e- x o 



<^(X)( X £(P))* 

- ^- L u {x))*xl u} (P) 



(C14) 



We kept the de Sitter modes in P-space because only 
|P| > k/A contribute to Eq. (69). Using Eq. (Cll), 
inserting the above expression in Eq. (69), and comparing 
the resulting expression with Eq. (C9), we get 



In Eq. (69), we need (the [/-mode contribution of) (2^ + 1)= / dP 1 dP 2 x R * {Pi)Xr( P 2) (C15a) 

,1 / tr T-» \ • 1 1 I v I ^ • _ • 1 1 _ l •* 



G" et (X,Pi) with \X\ X trans , since we are interested 



in the far away behavior of G^ c , and with Pj > v kA, be- 
cause the integrand vanishes for lower values of P. Since 
P U (X) -C Pi, the retarded character of Eq. (31) is auto- 
matically implemented, which means that 



x e - x ^- x ° 2 N UJ (P 1 ,P 2 ) 



GZ t (X,P 1 ) = (-i)G%(X,P 1 ) 



(Cll) 



The commutator G^(X,Pi), on the one hand, obeys 
Eq. (CI) in X, and on the other hand, behaves as in 



2c w = J dV 1 dV 2XR *(Pi)xT J *{V2) (C15b) 
x e~ x ^~ x ° > N U '(P 1 ,P 2 ) . 

These expressions are identical to those evaluated in 
de Sitter. Hence, n w and c w are respectively given by 
Eqs. (57a) and (57b). Therefore, to leading order in k/A, 
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and for an environment at zero temperature, and c u 
retain their standard relativistic expressions. 



-.Unruh 



This means that the state of the outgoing modes when 
they leave the central region III, and propagate freely, is 
the Unruh vacuum [34, 55, 56]. This can be explicitly 
checked from Eq. (C15) by re-expressing the out modes 
Xfl/i hi terms of the Unruh modes of Eq. (48). In this 
case, one finds that the mean number of Unruh quanta 



is given by, see Eq. (55), 



(2n 



Unruh 



l = 



dP 1 dP 2 {<t> bJ (P 1 ))*^{P 2 ) 



_ ( C16 ) 



x e-J-o -A) N W (P U P 2 ) 
= 1 + 0(k/A) . 

In other words, the role of the double integrals in 
Eq. (C15) and Eq. (C16), whose integrand explicitly de- 
pends on the actual "trans-Planckian" physics governed 
by A, f(P), r(P), is to implement the Unruh vacuum in 
dissipative theories. 
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